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The solution of non-linear evolution equations for dense nuclear gluon density has been suggested 
as one of the relevant mechanisms of pA and AA collisions at collider energies. Here we study a simple 
parameterization for the unintegrated gluon distribution using the knowledge of asymptotic solutions 
of the Balitsky-Kovchegov equation, describing high-energy QCD in the presence of saturation 
effects. A satisfactory description of nuclear shadowing at small- a; is obtained and it allows us to 
understand the qualitative behavior shown by data. 
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I. INTRODUCTION 



The knowledge of the QCD dynamics at high energies 
is essential in understanding the hadronic interactions 
studied at current (HERA, Tevatron) and future (LHC) 
accelerators. In the physics of saturation [1, „.2, 3, 4, 5., Q 
the relevant observable is the forward scattering ampli- 
tude of a quark- antiquark dipole off a target. It enters 
several processes at high energy like deep inelastic scat- 
tering (DIS), photoproduction and hadron-hadron scat- 
tering. Its quantum evolution has recently been the ob- 
ject of many studies J^[l,[l,[l3|- The Balitsky-Kovchegov 
(BK) equation [ill, valid in the large Nc limit and in 
the mean field approximation provides a tool to study 
the rapidity behavior of the gluon distribution in a large 
momentum range including the saturation region where 
the effects of gluon recombination, taken into account 
in the non-linear term, become important. The main 
features of anal ytic al solutions of BK have been used 
to parameterize |14| the unintegrated gluon distribution, 
tp{k,Y), which is the Fourier-transform of dipole-target 
amplitude. The knowledge of the corresponding nuclear 
unintegrated gluon distribution is indeed crucial. It is the 
basic object required to calculate physical observables in 
the A:j_-factorized picture employed by the gluon satura- 
tion models. In principle, ip{k, Y) possesses a Bjorken-a; 
dependence determined (here, Y = log(l/a;)) by non- 
linear evolution equations of the CGC theory [ij] and 
its dependence on transverse momentum k is fixed by a 
characteristic saturation momentum, Qs(Y). In most of 
CGC approaches, the gluon distribution is suppressed be- 
low the saturation scale tpA oc log{Ql/k^) [l^ compared 
to the perturbative form (pA oc k^^ . The physics of dense 
partonic systems and their non-linear perturbative evolu- 
tion to higher energy has motivated several attempts at 
understanding bulk properties of ultra-relativistic heavy 
ion collisions such as the multiplicity, rapidity distribu- 
tion e centrality dependence of particle production [l^ . 
Notice that saturation effects has been claimed to ac- 
count for the suppression of the high-py hadronic spec- 
tra in gold-gold and deuterium-gold collisions at RHIC. 



For the evolved gluon distribution, the yield of produced 
gluons in pA and AA collisions at central rapidity can be 
calculated in the factorized way [l|. 
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where k' = [k — pt)- Under the assumption of local 
parton-hadron duality, the multiplicity in AA collisions 
at central rapidity rises proportional to the nuclear sat- 
uration scale, dNAA/dr] oc Ql Aiv ^ 0) ^or col- 
lision energies reaching to 5.5 TeV, the upcoming pro- 
gram in lead-lead collisions at the CERN Large Hadron 
Collider (LHC), the confirmation of the conclusions ob- 
tained from RHIC is expected to take place. It is also 
expected the discrimination between the different physi- 
cal mechanisms proposed to explain particle production 
in high energy nuclear reactions. As noticed in Ref. [l6[, 
purely empirical parameterizations of multiplicity data 
of a large variety of colliding systems allow a logarithmic 
dependence on collision energy. However, RHIC data by 
themselves do not differentiate between this and other 
functional forms like power-laws, negating any possibility 
to usefully constraint the expectations for LHC energies 
without further theoretical guidance. 

In this work we propose a parameterization for 
the nuclear gluon distribution based on the numeri- 
cal/analytical solutions of the BK evolution equation. 
Starting with a a proton target, we present a parame- 
terization which encodes the main features of analytical 
solution of BK equation as geometric scaling property 
and inclusion of sub-leading terms. It describes all data 
on DIS in the region of a; < 0.01 in a wide interval of pho- 
ton virtualities. Relying on geometric scaling arguments, 
we propose a simple model for the nuclear unintegrated 
gluon distribution. Such a model is then compared to 
data on nuclear ratio and its relation to saturation mod- 
els is investigated. 
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II. A MODEL FOR THE UNINTEGRATED 
GLUON FUNCTION 

Lets start by considering the coUision between a virtual 
photon and a proton at high energy. In a frame where 
the photon travels fast, one can consider that it fluctu- 
ates into a qq dipole 17]. The lifetime of this dipole being 
much longer than the time of interaction with the pro- 
ton, one can write the cross section as a product of the 
wavefunction for a photon to go into a dipole times the 
dipole-proton cross section, which leads to usual formula 



dip o 

m 



(fr 



dz |vI/i.,i(r,z;Q^)|^a]i;(r,F), 

(1) 



where o'^ip^(f', ^) is the dipole-proton cross section. 
The transverse and longitudinal photon wavefuctions in 
this expression are computable in perturbative QED. 
The proton structure function F2 can be obtained 
from the 7*p cross section through the relation F2 — 
Q'^ / {■Ke^)[aT + <Jl]- 

If one treats the proton as a homogeneous disk of ra- 
dius Rp, the dipole-proton cross section in Eq. ([1]) is 
usually taken to be proportional to the dipole-proton for- 
ward scattering amplitude J\f{r, Y) through the relation 
<j2l^{r,Y) = 2nRlM{r,Y). Here, Y = log(l/x) is the 
rapidity variable. The BK evolution equation describes 
the high-energy evolution of the dipole-proton scattering 
amplitude Af, where the asymptotic behavior of its solu- 
tions is naturally expressed in momentum space. There- 
fore, the photo-absorption cross section can be written 
in terms of (p{k, Y), the Fourier transform of JV{r, Y): 
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(2) 



Accordingly, the proton structure function can be ex- 
pressed in momentum space: 



em Jo 
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where the wavefunction is now expressed in momentum 
space as well (see Ref. for details). The scattering 
amplitude ip{k,Y), the forward scattering amplitude in 
momentum space, can be obtained as a numerical result 
of the BK evolution equation. In case of neglecting its 
dependence on the impact parameter, the BK equation 
can be expressed as 

9y^(fc, Y) = axi-dL) ifiik, Y)~a ^\k, Y), (3) 

where a — asNc/ir and x{l) = 2'!/'(l) ~ "0(7) ~ "0(1 ~7) is 
the characteristic function of the Balitsky-Fadin-Kuraev- 
Lipatov (BFKL) kernel [18] and we used L = log(A:^/fcg) 
with fco some fixed soft scale. 



In Ref. [1J|, the properties of analytical solutions of 
BK equation were used to build an analytical expression 
for (p{k, Y). For instance, the BK equation reduces to the 
Fishcr-Kolmogorov-Petrovsky-Piscounov (F-KPP) equa- 
tion 19] (after changing of variables) when its kernel is 
approximated in the saddle point approximation i.e. to 
second order in the derivative Ol, the so-called diffusive 
approximation. The F-KPP equation is a well-known 
equation in non-equilibrium statistical physics, whose dy- 
namics is used to describe many reaction-diffusion sys- 
tems in the mean-field approximation. This equation has 
been extensively studied recently and, in particular, it is 
known to admit traveling waves as asymptotic solutions 
[lol ]. At asymptotic rapidities, the amplitude tp{k, Y), in- 
stead of depending separately on k and Y, depends only 
on the scaling variable Tp = k^/Q^, where we have intro- 
duced the saturation scale Ql{Y) oc exp{vcY), measuring 
the position of the wavefront. As discussed in Ref. 
a more detailed calculation allows also for the extraction 
of two additional sub-leading corrections, resulting into 
the following expression for the tail of the scattering am- 
plitude: 



ipp {k, Y) « (Tp) log (Tp) exp 



log' (Tp) 



2ax"ilc)Y 



, (4) 



where x" denotes the second derivative of the BFKL ker- 
nel with respect to 7 and where the saturation scale con- 
tains sub-leading terms. The critical parameters 7c and 
Vc are obtained from the knowledge of the BFKL kernel 
alone and correspond to the selection of the slowest pos- 
sible wave, i.e. Vc = ax'ilc)- For the LO BFKL kernel, 
one finds jc — 0.6275 and Vc = 4.88a. 

Concerning the property of geometric scaling [20j . 
when one moves along the saturation line, the behav- 
ior of the scattering amplitudes remains unchanged. In 
addition, the sub-leading corrections in ^ also play an 
important role as the last term introduces an explicit 
dependence on the rapidity Y and hence violates geo- 
metric scaling. However, this term in negligible when 

2s'x' i'y')Y ^ ^' ^^^^ means that geometric scaling is 



obtained for log (rp) < y/2x"{jc)aY, i.e., in a win- 
dow which extends Vy above the saturation scale. It 
is a remarkable property that, at high-energy, the con- 
sequences of saturation are observed arbitrarily far in 
the tail, where (p(k,Y) is much smaller than 1. No- 
tice that Eq. (HJ only gives a description of the tail of 
the wavefront (fi{k,Y) <C 1 (fc ^ Qs). In the infrared 
domain, one can show that the amplitude behaves like 
(Pp{k ^ Qs) oc c — log {^/T^), with c being a constant. 

In Ref. [l3|, an analytic interpolation between both 
asymptotic behavior of the amplitude was proposed. The 
final expression for the unintegrated gluon function is 
given by: 



ipp I 
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log 



[1 - exp (-T, 
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(5) 
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This equation is our master formula for the further ex- 
trapolation to the nuclear case. In expression above, T^n 
is an expression which is monotonically decreasing with 
L and which reproduces (up to the logarithmic factor), 
the amplitude for geometric scaling (|4]): 



TdiiiTp,Y) = exp 



-7c log (rp) 



log2(l+Tp)-log2(2) 



2ax"{lc)Y 



(6) 



with = 



2 ^aVaY 



The parameters of model have 



been determined from a fit of all the last HERA mea- 
surements of the proton structure function from HI [2]| . 
ZEUS [l^l, with the analysis being restricted to the fol- 
lowing kinematic range x < 0.01 and 0.045 < < 
150 GeV^. In our numerical analysis we use the fit in- 
cluding charm [ll| (m, = 0.14 GeV and mc = 1.3 GeV), 
where = 0.807, x" = 2.96, = 3.917 x IQ--"' GcV^ 
and Rp = 4.142 GeV. 

Having presented the unintegrated giuon distribution 
for the nucleon case, in what follows we investigate its 
extension for nuclear targets relying on geometric scaling 
properties. Based on geometric scaling arguments and on 
the main features of numerical solutions of BK equation, 
in Ref. [2^ it is assumed that both the energy and the 
nuclear size (or centrality) dependence on the scattering 
amplitude N(r, Y; b) can be encoded in the saturation 
scale Qs,hiY,b) for any hadron h (proton or nucleus). 
Then the cross section in Eq. ([T]) can be written as 



Y) - 7:Rl 



dz |^'T,L(r,z;Q2)| 



X 2 / cPhUhirQ.x, h)- 



(7) 



As the wavefunction is proportional to times a 
function of r^Q^, the cross section is only a function 
of T/i = Q'^/Qly^. This geometric scaling was found to 
describe all lepton-proton data with x < 0.01. In or- 
der to compare protons and different nuclei, in Ref. [l^ 
the following assumption about the impact parameter 
dependence was considered: it can be scaled by the nu- 
clear radius of the hadronic target h,b = 6/-y/7rR^, with 
Ra = (1.12^1/3 „ o.86A~i/3) fjjj^ Tj^g^^ the condition 
for geometric scaling in lepton nucleus data is 




fnR\ 



\TrRl 



(8) 



For the dependence on atomic number, the saturation 
scale in the nucleus grows with the ratio of the transverse 
partons densities to some power A, which was taken as 
a free parameter: 



' AttRI' 



(9) 



where QJY) is the saturation scale for a proton target. 
In Ref. it was found A = 1/5, with 5 = 0.79 ± 0.02, 



1.2 

1.1 

Q 1.0 
O 

O 0.9 
< 0.8 
0.7 




1.2 
1.1 - 

f 1-0 
O 

O 0.9 
< 0.8 

DC 

0.7 




FIG. 1: The nuclear ratio R(A/B) = BF^/AF2 as a function 
of Bjorken-a; for several nuclear targets. Solid line represents 
our results (BM model) and dot-dashed curves the result of 
ASW model. 



which translate into a growth of the nuclear saturation 
scale faster than A^^^ for large nuclei. 

Now, we are ready to construct a model for the unin- 
tegrated nuclear gluon distribution. Using the geometric 
scaling arguments, we can rewrite Eq. ([5]), by replacing 
Rp ^ Ra and Q,(y) ^ QsAY)- 

^A{k,Y)^ log(^I^^+l {l-cxp[-rdii(rA,r)]}, 

(10) 

where we have defined ta = {R\/ARl)'^Tp. We think 
this procedure is suitable once it has been shown in 
Ref. [2^ that the data for nuclear structure functions 
lie on the region where ta < 10. Therefore, the condi- 
tion log {ta) < y/2x"i-/c)aY ~ 1.22^1og(l/a;) is roughly 
valid and geometric scaling violation should be small. In 
next section, we investigate the numerical calculation of 
the nuclear structure function using Eq. (|10p and com- 
pare them to the experimental measurements of the nu- 
clear ratios. 



III. DISCUSSIONS AND CONCLUSIONS 

Let us now compare our model for the nuclear un- 
integrated gluon distribution to the experimental data 
on nuclear ratios. The shadowing in nuclei is usu- 
ally studied through the ratios of cross sections per nu- 
cleon for different nuclei. In Fig. [1] the results of the 
model are compared with experimental data from Refs. 
dil. In these figures, the quantity RAT10{A/B) = 
[BF2Aix,Q'^)/AF2Bix,Q'^)] is presented. We have se- 
lected data where x < 0.01 and only statistic errors are 
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shown. The sohd hnes (here labeled as BM model) rep- 
resents the numerical results using A = (0.79 ± 0.02)^^. 
This means that we are using the conservative value of 
A which reproduces the same growth for the nuclear 
saturation scale as in Ref. [1^. It is verified that we 
fairly reproduce the behavior on x for the nuclear ra- 
tio, predicting a larger nuclear shadowing for the carbon 
case. The model is able to describe data even in the re- 
gion X < 10^^, taking part of the anti-shadowing region. 
Considering both the absence of any free parameter in 
the nuclear unintegrated gluon distribution we find the 
agreement quite reasonable. For sake of comparison, we 
show the numerical results using the Armesto-Salgado- 
Wiedemann (ASW) model [1^ (dashed lines) for the nu- 
clear cross section. It is verified a saturation of nuclear 
ration for x > 10~^ and the model underestimate the 
data for Lead target. For consistency, we have checked 
the results for NMC data {x > 0.0125) and a smaller 
A values seems to be preferred. This fact has been dis- 
cussed in Ref. [25| , where the nuclear enhancement in the 
saturation scale has been addressed in distinct saturation 
model (with and without geometric scaling violation). As 
pointed out in that study, a precise extraction on the A 
dependence of Qs will play an important role in distin- 
guishing between "classical" and "quantum" evolution in 
the CGC. Similar discussion can be found in Ref. ^26]. 

It is timely discuss the qualitative aspects of the nu- 
clear ratios using the present parameterization for the un- 
integrated gluon distribution. Assuming geometric scal- 
ing in the region covered by the experimental data, we 
can investigate the asymptotic behaviors of the nuclear 
ratio. At intermediate x, around x = 0.1 the unin- 
tegrated gluon distribution has the power-like behavior 
(modulo logarithmic corrections) ipA oc {Ql a/^'^V" ^^"^ 
(fip oc {Ql/k'^)'^'' (in fact, for the proton case we have ad- 
ditional scaling violations). This leads to the following 
result for the nuclear ratio, taking into account geometric 
scaling arguments: 



ARl 



l-7cA 



,(11) 



where the asymptotic behavior seems to be dependent 
on the product 7cA. Notice that in our case 7c = 0.63 
against 7c — 0.75 for the ASW parameterization. There- 
fore, the expression above help us to understand the de- 
viations between ASW model and the current model. We 
have checked that Eq. (fTTj) exactly gives the plateau in 
R{A/B) seen in the ASW results in the region x > 0.01. 
The enhancement observed in our case is a consequence 
of the anomalous dimension to be dependent on x and 
in that kinematic region, where scaling violations start to 
take place. Therefore, the discussion is more involved in 
that region, which resembles the discussion on the Cronin 
peak in A A collisions. 

On the other hand, for very small- a; (and low Q^), 
where k < Qs{Y), the unintegrated gluon distribution 



takes a logarithmic behavior. Thus, one has (pA oc 
^ogliQlA/k^V"] and ifp oc logHQl/k^)-^"]. This leads 
to the following result for the nuclear ratio. 



AFi °^ i?2 



"l0g(Q?,A)' 


R\ 


_ log(Q?) _ 


ARl 



iog(§;) 



(12) 



which is clearly dependent on the saturation scale and 
on photon virtuality. In particular, the ratio seems to be 
closely independent of anomalous dimension 7c. 

As a final discussion, we discuss the relation of present 
model with recent determinations of the nuclear unin- 
tegrated gluon distribution. Notice that the usual nota- 
tion for the unintegrated gluon density, /(fc, Y) is related 
to function (p(k, Y) by the following relation (see for in- 
stance Ref. [27] for a derivation): 



f(k,Y) 



Nc 



i-fc2_l) ep{k,Y). (13) 



There is an comprehensive research program in calculat- 
ing the numerical solution for the nuclear gluon distribu- 
tion from the BK equation including NLO effects such as 
the running coupling corrections [l6|, [2^ . From the phe- 
nomcnological point of view, our results are similar to 
the numerical solution as on average the energy behavior 
for the saturation scale is consistent with A = 0.2 and 
the main features of solutions are present. Recently, a 
new modified BFKL equation incorporating the shadow- 
ing and anti-shadowing corrections of the gluon recombi- 
nation to the unintegrated gluon function has been pro- 
posed through the RSYZ equation [2^. The numerical 
results from such an equation is not identical to the BK 
results in the transverse momentum space but is similar 
to a mean field result. Of course, our model is less general 
than gluon function obtained from RSYZ equation as we 
focus on the shadowing region. Finally, in Ref. [30] the 
authors have tested several parameterizations of uninte- 
grated gluon function in nucleus inspired on asymptotic 
solution of BK. Most of them fit our expression for the 
ip{k,Y), which gave satisfactory description of the data 
for ratio central/peripheral at RHIC. 
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